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ABSTRACT 



Based on high quahty spectral data (spectral resolution R>60000) within the wave- 
length range of 3550-5000 A we determined main parameters (effective temperatm'e, 
surface gravity, microturbulent velocity, and chemical element abundances including 
heavy metals from Sr to Dy) for 14 metal-deficient G-K stars with large proper mo- 
tions. The stars we studied have a wide range of metallicity: [Fe/H] = —0.3 -. — 2.9. 
Abundances of Mg, Al, Sr and Ba were calculated with non-LTE line-formation effects 
accounted for. Abundances both of the radioactive element Th and r-process element Eu 
were determined using synthetic spectrum calculations. We selected stars that belong 
to different galactic populations according to the kinematical criterion and parameters 
determined by us. We found that the studied stars with large proper motions refer to 
different components of the Galaxy: thin, thick disks and halo. The chemical composi- 
tion of the star BD+80° 245 located far from the galactic plane agrees with its belonging 
to the accreted halo. For the giant HD 115444 we obtained [Fe/H] = —2.91, underabun- 
dance of Mn, overabundance of heavy metals from Ba to Dy, and, especially high excess 
of the r-process element Europium: [Eu/Fe] = -1-1.26. Contrary to its chemical compo- 
sition typical for halo stars, HD 115444 belongs to the disc population according to its 
kinematic parameters. 

Subject headings: high resolution spectroscopy - ground-based UV and blue region - 
metal-poor stars - chemical composition 



1. Introduction 

Spectral and kinematical parameters of unevolved stars that belong to different galactic popu- 
lations are the main sources of information about chemical and kinematical evolution of the Galaxy. 
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We see that spectroscopy of stars of different generations in the Galaxy is a key program for the 
largest telescopes in the world. During last decades scientists become more and more interested 
in spectroscopic studies of the oldest stellar populations, since chemical evolution of the Galaxy is 
imprinted in the chemical composition of metal-poor stars of subsequent generations. To recon- 
struct the chronology of the chemical and kinematic evolution of the Galaxy we need an enormous 
number of high quality spectral data and accurate proper motion and parallax determination for 
metal-deficient stars. 

Unevolved F-G stars (dwarfs and subdwarfs) are the most effective probes of the chemical 
composition and measurements of radial velocities. Their spectra are abundant in narrow and 
low-blended absorptions as seen at a high spectral resolution. In order to enlarge the sample of 
unblended lines, we have to record a broad spectral region. For extremly metal-deficient stars that 
have weak lines, the spectra have to be taken in the blue and ultraviolet regions where there are 
more lines than in the visible range. Considering the facts mentioned above, we performed high 
resolution spectroscopy of a sample of F-G metal-poor stars at the 6-meter telescope of the Special 
Astrophysical Observatory of the Russian Academy of Sciences. All our spectra were taken within 
a broad spectral interval from 3550 to 5000 A. Below we present some recent results. 



2. Spectral observations and spectra reduction 



We have selected a sample of high- velocity stars from the high proper motion stars review (jCarnev et al 



1994 ). HD/BD numbers o f 14 stars are listed i n Table[TJ The spectral data were taken using the 
echelle spectrograph NES (jPanchuk et al.ll2009l ) permanently mounted in a Nasmyth focus of the 
6-meter telescope of the Special Astrophysical Observatory. NES provides a spectral resolving 
power of R>60000 in the spectral range 3000-10000 A. NES is capable to observe in the UV due to 
the camera made of fused silica. Now NES works in combination with a 2048x2048 pixel CCD that 
has a high sensitivity in the ground-based UV spectral range, for wavelengths longer of 3000 A. We 
notice that spectral range of the spectrog raph NES is being crosse d with the range (1740-3100 A) of 
the high resolution spectrograph UVES (jKappelmann et al.ll2006l ) which is pr ojected for the World 
Space Observatory-UV (more details see in the paper by Shustov et al. l2009l ) . 



The 2D echelle-spectra were reduced (applying standard procedures of bias subtraction, scat- 
tered light and cosmic ray trace removal, and order extraction) using the context ECHELLE of 
the system MIDAS. The conte xt ECHELLE was moderniz ed to process spectra obtained in com- 
bination with an image sheer (jYushkin &: Klochkoval |2004| ) . The signal-to-noise ratio for all the 
spectra discussed in this paper is higher than 200. Combined with the spectral resolving power, 
that allowed us not only to detect rather weak lines but also to study their profiles. 

Kinematic characteristics (pr oper motions and parallaxes) of th e stars we study are presented 
in Table 1 in our previous paper (jKlochkova &: Tavolganskayall20ld ). In this paper, based on the 
same spectra taken with the 6-meter telescope, accurate radial velocities of 15 metal-deficient stars 
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were published (jKlochkova &: Tavolganskavall201C ). The list of spectral lin es (more than 8000) was 



made using the VALD database (jPiskunov et al.l Il995l : iKupka et alj Il999l ) . Nearly 860 unblended 
features from this list were selected for Vr measurements. The standard deviation of the measured 
velocity does not exceed a < 0.9km/s for stars with the metallicity [Fe/H]> —1, and a < l.lkm/s 
for stars with [Fe/H]< —1. 



3. Model parameters and elemental abundances determination 



The effective temperature T(,ff was det ermined usin g the Stromgren uvbyl3 indices and cali- 
bration based on the infrared flux method (jAlonsdll996l ). Metallicity values needed for the first 
iteration of Tf,ff determination were used from publications, but in the following iterations our 
spectroscopic [Fe/H] values were used. In this procedure from 100 to 230 Fe i lines were taken 
into account for different stars. Additionally, to control the Tg// values we used the generally 
adopted spectroscopic way of Tg// determination, forcing independence of Ige(-Fe) on the lower 
level excitation potential. 



Surface gravity was calculated using known relations: 
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where M^oi = V + BC + 51g7r + 5, - mass of the star, M^oi is the bolometric luminosity, V is 
the visual magnitude, BC denotes bolometric correction, it is the stellar parallax. 



The Hipparcos parallaxes (jPerrvman et al .1119971 ) were used for th e calculations. Bolometric 
corrections were calculated using a calibration formula of Flower (|l99a ) . 



The microturbulent velocity £,t was determined forcing the independence of the neutral iron 
abundance on the equivalent width Wa of the line. Atmospheric parameters adopted for 14 stars are 
presented in Table[TJ The WIDTH9 code and Kurucz's grid of atmospheric models (|Kurucall993l ) 
were used to calculate chemical abundances. To obtain the normalized logarithmic abundance ratio 
of elements to iron [Elem/Fe] = lo g(A^rapm/-/VFR) — ^og{NEiem/^Fe)Q we used solar abundances 



from VALD (Piskunov et al. 



1995 



presented by Asplund et al. (120091). Spectral lines f or the calculation were selected using data 



Kupka et al.l Il999l ) . The atomic constants of lines needed for 



abundance calculati ons w ere taken from the same database. The list of lines was published earlier 
by Klochkova etal. (j2006l ) and it is available at http://w ww. chjaa. org/2006- 6- 5.htm\ 



The use of echelle spectra observed within a broad short-wavelength spectral range allowed the 
model atmosphere method to be applied to determine abundances of 20 chemical elements including 
heavy metals from Sr to Dy presented in TableO Abundances of the iron group elements (Cr, Mn, 
Co, Ni, Zn), which origin in the equilibrium Si-burning, show very small scattering relative to the 
Fe abundance. Especially it concerns zinc, whose average relative abundance is [Zn/Fe] =0.06 for 
stars with [Fe/H]< -1.4. 
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Fig. 1. — The Mgl A 4703 A line profile (top) observed in the spectrum of HD 148816 (dotted line). 
A synthetic profile is shown by a thick line 



Abundances of Mg, Al, Sr and Ba are calculated with a non-LTE approximati on (see Table[2|) . 
The non-LTE abundance elements was determi ned w ith a modified MULTI code ()Carlssonl ll986). 
Modifications are described by K orotin et al. (1 19991 ). Since we use Kurucz's atmosphere models 
calculated with ATLAS 9 (Kurucz 119931 . |2005| ) , the necessary background opacities for MULTI are 
also from ATLAS 9. In the modified version, the mean intensities that are used to obtain the 
radiative photoionization rates are calculated for a set of frequencies at each atmospheric layer, 
and then they are stored in a separate block, where they can be interpolated. 

The atomic model of magnesium used in this work is essentially the same as that described by 
Mishenina et al. (120041 ). This model consists of 84 levels of Mgl, 12 levels of Mgll and the ground 
level of Mglll. Within the described system, transitions between first 59 levels of Mgl and ground 
level of Mgll have been considered. The detailed structure of the multiplets was ign ored and each 
LS multiplet was considered as a single term (more details in lMishenina et al.ll2004l ). 

For the non-LTE calculations we adopted an aluminium atomic model that consists of 78 
levels of All and 13 levels of AlII. Oscillator strengths and photoionization cross sections were 
taken from the TOPbase ( http://vizier.ustrasbg.fr/topbase/topbase.html) . More details see in the 
paper by Andrievsky et al. (120081 ). 

Since the resonance Sr II and Al I lines in the UV band of the spectrum are somewhat blended 
with other metallic lines, to compare non~LTE line profiles with the observed spectrum one needs 
to use a combin ation of a non- LTE and LTE synthetic spectrum. This was done with the updated 
code SYNTHV (jTsvmballll996l ) which is designed for synthetic spectrum calculations within LTE. 
With this program we calculated synthetic spectra for selected regions comprising the SrII and 
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Fig. 2. — The same as in Fig. 1, but for the Hne Ball A 4554 A 



All hnes of interest taking into account all t he lines from each region listed in the VALD data- 



base (Piskunov et al 



1995 



Kupka et al.lll999l ). Then we included corresponding 6-f actors (factors 



of deviation from LTE level po pulations) calcu lated in MULTI separately for the strontium and 
aluminium lines into SYNTHV (|Tsymba]|ll996l ). where they were used for calculations of the non- 
LTE line source function. 

Our strontium atomic model consists of 44 levels of Sr II with n < 13 and Z < 6, as well as the 
ground level of SrIII. Since the ionization potential of the neutral strontium is only 5.7 eV, even 
in atmospheres of cool stars it exists only as of SrII. Therefore, 24 levels of Sri were included in 
the model only to conserve the particle number. Radiative photoionization rates for the s, p and d 
levels are based on th e phot oionization cross sections calculated with the quantum defect method 
(see Andrievsky et al. I2OIOI for details). 



Our barium model contains 31 levels of Bal, 101 levels of Ball with n < 50, and the ground 
level of Ba III ion. We considered 91 bound-bound transitions in detail. A cause of some uncertainty 
in the non-LTE analysis of the barium spectrum is the scarce information about the photoionization 
cross-sections f or different levels. We used results obtai ned u sing the scaled Thomas-Fermi method 
(|Hofsaesall979l ) (for more details see Andrievsky et al. |2009| ). 



To illustrate these calculations, we compare the Mg 1 4702 A profile observed in the HD 148816 
spectrum with a theoretical one in Fig. 1. A similar comparison for the lines Ball A 4554 A, 
Eu II A 4129 A and Th II A 4019 A is presented in Fig. 2-4. 

Aiming to estimate cosmochronological age, the abundance ratio of the radioactive element 
Th to that of the r-process element Eu is determined for the selected stars (see Tabled]). Being 
independent on the stellar evolution models, the cosmochronological method of the age estimation 
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is sensitive both to reliability of stellar nucleosynthesis calculations and to accuracy of heavy metal 
abundances derived from spectroscopic observations. We determined Eu and Th abundances by 
comparison of the observed and synthetic spectra in the region of the 4129 A and 4019 A lines for 
Eu and Th, respectively. Synth etic sp ectra were calculated using a current version of the STARSP 
package developed by Tsymbal ^9M)- The Th4019 A line is heavily blended (Fel, Nil, Col, VII, 
^^CH), therefore the error of the Th abundance determination is over 0.1 dex. Taking into account 
these errors in the abundances of Eu and Th, we (as well as other authors) obtained too large age 
errors, more than 4 Gyr. For this reason we do not present our age estimations. 



Discussion of results 



The stars were related to subpopulations of the Galaxy on the basis of the components of stellar 
spatial velocities. Then this stratification was refined using elements of Galactic orbits of the stars. 
To do this, we computed the orbital elements by simulating 30 revolutions of the star around the 
Gala ctic center, based on the model of the Galaxy containing a disk, bulge, and extended massive 
halo (|Allenlll99ll ). The Galactocentric distance of the Sun was assumed to be 8.5 kpc, the rotational 
velocity of the Galaxy at the solar Galactocentric distance was 220km/s, and the velocity of the Sun 
with respect to the local standard of rest is (U,V, W)© = (— 11, 14, 7.5) km/s (IRatnatunga et al 



19891 ). The stars of the thin and thick disc subsystems were separated according to a method that 
calculates the p robability of the star membership in each subsystem using components of their 
spatial velocity (jKoval et al.l l2006l ) . Stars of the thick disk and proper halo were separated using 
total spatial velocity relative to the local standard of rest - Visr = 175km/c (jMarsakov. Borkova 
20051 ). The accreted halo stars (i.e. stars of extragalactic origin) were identified using criterion 
- Visr > 240km/c - these stars have high orbital energies and most of them move on retrograde 



orbits (jMarsakov. Borkoval l2006l ) . 



In this work, the parameters (including the metallicity) and elemental abundances were deter- 
mined for 10 stars. We also calculated the Mg, Al, Sr, Ba abundances in a non-LT E approximation 
for 4 stars using parameters that were determined earlier iKlochkova et al.l (|2006l ) . As follows from 
Tabledl the stars have metallicities within a wide interval: [Fe/H] = —0.3 -. — 2.9. The relative 
abundances [Elem/Fe] presented in Tables[3l and [2 ] correspond to the expected values according to 



common behaviour from papers (jMcWilliam 



1997 



Roedereii l2009l ) . This means that, for example, 
the a-process elements (Mg, Si, Ca, Ti) are overabundant for metal-poor stars at [Fe/II]< —1.0, 
and the value [a/Fe] is almost constant with a small dispersion: [a/Fe]= +0.3 + 0.4. For stars at 
[Fe/II]> —1.0 abundances of these elements decrease. The star BD+80°245 is an exception. 

Our resul ts for [Mg/Fe] are very well correspondent to the [Mg/Fe]-behaviour derived by 
Gehren et al. (|2006l ) for an extensive stellar sample. 



Taking into account kinematical parameters and chemical composition of stars studied, we rec- 
ognized their belonging to various galactic populations. The most noteworthy result is a peculiar 
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Fig. 3. — The same as in Fig. 1, but for the Hne EuII A4129 A in the spectrum of HD 22879 
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Fig. 4. — The same as in Fig. 1, but for the hne Th II A 4019 A in the spectrum of HD 5256. Synthetic 
spectra (from top to bottom) are calculated with N(Th) = —12.5, —12.0, —11.5 

chemical abundance pattern in the atmosphere of BD+80°245 whose metallicity is [Fe/H]= —1.71. 
The star is allocated with both lowered (for its metallicity) abundances of Mg, Al, Ca and remark- 
able deficiency of the s-procees element Ba: [Ba/Fe]= —1.46. As follows from Tabled BD+80°245 
is also very distant from the Galactic plane Zmax = 23kpc. Summarizing our results (metallicity, 
peculiar chemical composition, and kinematical parameters), we classify this star as belonging to 
the accreted halo. It should be noted that this conclusion is in agreement with earlier results of 
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Carney et alj (|l997l ) who revealed similar chemical peculiarities for BD+80°245 and suspected that 
it is an object originated not in our Galaxy. Now it is known that the galactic halo is not a single 
homogeneous population, but at least two subpopulations, outer and inner, (or initial and accreted) 
(|Venn et alJbood : ICarollo et allboiol ). 



It should be noted that other two stars, BD+71°31 and HD 148816 which kinematically belong 
to the accreted halo, do not have chemical peculiarities similar to those obtained for BD+80°245. 

There is another paradoxical object - a giant HD 115444 with the low metallicity, [Fe/H] = — 2.91. 
For this star we obtained a pecul iar ch emical abundance pattern which is in agreement with pub- 
lished earlier by Westii i et al. (|2000l') : overd eficiency of Mn and high overabundance of heavy 
metals from Ba to Dy (iKlochkova et al.ll2006l ). We noted especially high excess of the r-process 
element Europium: [Eu/Fe] = +1.26. On the one hand, we see a star with low metallicity and large 
overabundance of heavy metals. In total, the abundances of chemical elements are typical for halo 
stars. On the other hand, kinematically HD 115444 belongs to the disk population. It is evident 
that its origin and history do not correspond to each other. 

More detailed results concerning chemical abundances, classificati on by population, and ele- 
ments of cosmochronology will be presented in the forthcoming paper (IKlochkova et ahlbnill ). 



Summary 

Based on spectroscopic data (R>60000 within the wavelength range 3550-5100 A) taken with 
the spectrograph NES of the 6~meter telescope, we determined parameters (effective temperature, 
surface gravity, microturbulent velocity, metallicity, and accurate radial velocity) of a sample of 
14 metal-deficient G-K stars with large proper motions. The use of echelle spectra allowed the 
model atmosphere method to be applied for determination of abundances of 20 chemical elements 
including heavy metals from Sr to Dy. Abundances of Mg, Al, Sr and Ba are calculated within a 
non-LTE approach. 

Taking into account kinematical parameters and chemical composition of the investigated 
stars, we determined galactic populations they belong to. The chemical composition of the halo 
star BD+80° 245 that is located far from the galactic plane agrees with its belonging to the accreted 
halo. In addition to the overdeficiency of a-process elements, the star has the underabundance of 
the s-process element barium: [Ba/Fe] = —1.46. Peculiarities of the chemical composition derived 
for another metal-deficient star HD 115444 do not correspond to its kinematic parameters. This 
star also has a non-standard origin. 

We conclude that the sample of high proper motion stars we studied is inhomogeneous - it 
includes objects within a wide metallicity interval [Fe/H] = —0.3 — 2.9 that belong to both disk 
(thin and thick) and halo (inner and accreted). 



-9- 



Acknowledgements 

V.P. and N.T. are much indebted both to the Russian Foundation for Basic Research for finan- 
cial support of this work (project 07-02-00247 a) and the Russian Federal program "Observational 
manifestations of evolution of chemical abundances of stars and Galaxy". T.M. and S.K. are grate- 
ful for the support of the Swiss national Science Fund (project SCOPES No. IZ73Z0-128180/1). 
V.M. is grateful for the financial support of the Ministry of education and science of the Russian 
Federation (the project P. 685). 



-10- 



REFERENCES 

Allen, C, Santillan, A. 1991, Rev. Mex. A&A, 22 255 

Alonso, A., Arribas, S., Martinez-Roger, C. 1996, A&A, 313 873 

Andrievsky, S.M., Spite, M., Korotin, S.A., et al., 2008, A&A, 481 481 

Andrievsky, S.M., Spite, M., Korotin, S.A., et al., 2009, A&A, 494 1083 

Andrievsky S.M., Spite F., Korotin S.A., et al., 2010, A&A, 509 88 

Asplund, M., Grevesse, N., and Sauval, A.J., Scott, R 2009, ARA&A, 47 481 

Carney, B.W., Latham, D.W., Laird, J.B., et al., 1994, AJ, 107 2240 

Carney, B.W., Wright, J., Sneden, C, et al., 1997, AJ, 114 363 

Carlsson, M. 1986, Uppsala Obs. Rep., 33 

Carollo, D., Beers, T., Chiba, M., et al., 2010, ApJ, 712 692 

Flower, P. 1996, ApJ, 469 355 

Gehren, T., J.R. Shi, H.W. Zhang, et al., 2006, A&A, 451 1065 
Hofsaess, D. 1979, Atomic Data and Nuclear Data Tables, 24 285 
Kappelmann, N., Barnstedt, J., Gringel, W., et al. SPIE, 2006, 6266 25 

Klochkova, V.G., Zhao, G., Ermakov, S.V., Panchuk, V.E. Chinese Astron. Astrophys., 2006, 6 
579. 

Klochkova, V.G., Tavolganskaya, N.S. Astrophys. Bull., 2010, 65 18 

Klochkova, V.G., Mishenina, T.V., Panchuk, et al., Astrophys. Bull., 2011, 66 No.l, in press 
Korotin, S.A., Andrievsky, S.M., Luck, R.E., 1999, A&A, 351 168 
Koval', V.V., Marsakov, V.A., Borkova, T.V. 2006, Astron. Rep., 53 785 
Kupka, F., Piskunov, N.E., Ryabchikova, T.A., et al., 1999, A&AS, 138 119 
Kurucz, R. 1993, ATLAS 9 stellar atmospheres program and 2kms— 1. CD-ROM 13. 
Kurucz, R.L. 2005, Memorie della Societa Astronomica Italiana Supplementi, 8 14. 
Marsakov, V.A., Borkova, T.V. 2005, Astron. Lett., 31 515 
Marsakov, V.A., Borkova, T.V. 2006, Astron. Lett., 32 604 



- 11 - 



McWilliam, A., 1997, ARA&A, 35 503 

Mishenina T.V., Soubiran C, Kovtyukh V.V., Korotin S.A. 2004, A&A, 418 551 

Panchuk, V., Klochkova, V., Yushkin, M., Najdenov, 1. 2009, Journal of Optical Technology, 76 87 

Ferryman, M.A.C., Lindegren, L., Kovalevsky, J., et al., 1997, A&A, 323 L49 

Piskunov, N.E., Kupka, F., Ryabchikova, T.A., et al., 1995, A&AS, 112 525 

Ratnatunga, K.U., Balicall, J. N., and Casertano, S. 1989, ApJ, 339 106 

Roederer, I. 2009, AJ, 137 272 

Shustov, B., Sachkov, M., Gomez de Castro, A.I., et al. 2009, Astrophys. Space Sci., 320 187 
Tsymbal, V.V. 1996, Ap&SS, 108 198 

Venn, K.A., Irwin, M., Shetrone, M.D., et al., 2004 AJ, 128 1177 
Westin, J., Sneden, C, Gustafsson, B., Cowan, J. 2000, ApJ, 530 783 

Yushkin, M.V., Klochkova, V.G. 2004, Preprint of Special Astrophysical Observatory, No. 206 1 



This preprint was prepared with the AAS LM^^X macros v5.2. 



- 12 - 



Table 1: Atmospheric parameters and metallicity [Fe/H] of program stars. Binary stars are marked 
in the last column 
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Table 2: Non-LTE abundances of Mg, Al, Sr & Ba in the atmospheres of stars studied, A(Elem). 
Hydrogen abundance A(H)=12.0 
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Table 3: Relative abundances of chemical elements [Elem/Fe] 
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Table 4: Abundances loge(Eu) and loge(Th), relative abundances [Eu/H] and [Th/H], and ra- 
tios [Eu/Fe] and [T h/Fe]. Solar abundances loge(Eu) = 0.51 and e(Th) = 0.02 are taken from 
(|Asplund et alJbnnd ) 
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5256 
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148816 


215065 


log e(Eu) 


0.30 
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0.00 
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0.30 
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0.30 


[Eu/H] 
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loge(Th) 
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0.00 
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0.00 
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0.00 
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-0.02 






[Th/Fe] 


0.46 
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Table 5: Distance d, spatial velocity components U,V,W, Visr and parameters of the stellar orbits 
calculated: e - orbital eccentricity, Zmax - maximum distance from the Galactic plane, Rmax - 
apogalactic radius. Galactic subsystems to which the stars belong is noted in the last column: 1 
- thin disk, 2 - thick disk, 3 - initial halo, 4 - accreted halo, 1-2 - the star rather belongs to the 
thin disk than to the thick one; 3-2 - the star rather belongs to the initial halo than to thick disk 
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0.26 
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0.14 
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0.8 


13.2 


3- 


^2 


215065 


29.37±0.5 
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